A Measurement of the D + Lifetime 
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A high statistics measurement of the lifetime from the Fermilab fixed-target FOCUS 
photoproduction experiment is presented. We describe the analysis of the two decay modes, 
Df — * <^>(1020)7r + and Df — > K* (892)° K + , used for the measurement. The measured lifetime is 
507.4±5.5 (stat.)±5.1 (syst.) fs using 8961±105 D+ 0(1O2O)tt + and 4680±90 D+ -> K*(892)°K + 
decays. This is a significant improvement over the present world average. 

PACS numbers: 13.25.Ft, 14.40.Lb, 14.65.Dw 



Experimental measurements of charm lifetimes have 
given guidance towards a theoretical description of strong 
interactions at low energy scales 0, Q- For example, 
the ratio of the Df and D° lifetimes gives information 
on the weak annihilation contribution to the total decay 
0- LH EM @ • Not only is this important for improving our 
theoretical understanding of strong interactions, but the 
theoretical tools used to calculate lifetimes are the same 
or similar to those used in other areas, for example to 
extract |U cs |/|U c( i| and |U c b| in charm and bottom decays, 
or to calculate the 6-particle lifetimes 0, . 



The data used were collected by the FOCUS collabora- 
tion in the 1996-1997 fixed-target run at Fermi National 
Accelerator Laboratory. The FOCUS spectrometer is an 
upgrade of the spectrometer used in the E687 photopro- 
duction experiment Q . In the majority of the data used 
for this measurement (Run Period B) the vertex region 
consists of four BeO targets and 16 planes of silicon strip 
detectors (SSD). Two of the SSD planes were placed im- 
mediately downstream of the second target, and two im- 
mediately downstream of the fourth (most downstream) 
target. About 20% of the Df decays used in this anal- 
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ysis were from data taken earlier in the data-taking run 
(Run Period A) without the four SSD planes in the tar- 
get region 0|- Momentum analysis was done using five 
multiwire proportional chambers and two magnets with 
opposite polarities. Three multicell threshold Cercnkov 
counters were used for particle identification [Til ]. 
_The D+ -> 0[lO2O)vr+, (£(1020) -» K~ K+ and D+ -> 
K*(892)°K+, K*(892)° -> K~tt + resonant decay modes 
|22| are used because they have much better signal-to- 
background than inclusive Df — ■> i\T — _ftT + 7r + decays. 

All i\"~_ft' + 7r + candidates are tested to see if they form 
a vertex with a confidence level greater than 1%. The 
candidate Df momentum vector is then projected to 
search for a production vertex with one or more tracks. 
As many tracks as possible are included in the production 
vertex so long as the vertex confidence level is larger than 
1%. The production vertex is required to be within one of 
the four targets. To largely eliminate non-charm back- 
grounds, the separation L between the production and 
decay vertices is required to be larger than 6er l where ox 
is the calculated error on L. 

The (/>(1020)tt+ (K*(892)°K+) decay mode candidates 
are required to have K + K~ (K~ir + ) masses within two 
sigma of the nominal 0(1020) (iT(892)°) mass. Also 
the magnitude of cos#* must be larger than 0.3 (0.6) 
for the 0(1O2O)tt+ (K*(892)°K + ) decay mode candi- 
dates, where 8* is the angle between the K~ and the tt + 
(K+) intfie 0(1020) (K* (892)°) center-of-mass frame. 
For the K*(892)°K+ decay mode, the K+K~ invari- 
ant mass must not be within two sigma of the 0(1020) 
mass to ensure statistically independent samples in the 
two decay modes. Each track in the K~K + n + candi- 
date combination must also satisfy a minimal Cerenkov 
particle identification criteria. To further reduce the 



D+ 



K~ 



contamination by an additional factor 



of 15 in the D+ -> K*(892)°K+ candidate sample, the 
kaon with the same sign as the pion is required to have a 
kaon probability that is larger than its pion probability 
by a factor of 20. 

The K~K + tt + invariant mass plots for data are shown 
in Fig. ^ The fits shown are to a Gaussian signal and a 
quadratic background function which yields 8961 ± 105 
0(1O2O)tt and 4860 ± 90 K*{892)°K reconstructed D+ 
decays. The lifetime analysis uses K~K + tt + candidates 
in a signal region (SR) that is within ±2cr m of the fitted 
Df mass, where a m — 9.7 MeV/c 2 is the fitted Gaussian 
width of the Dj peak. Candidates within two symmetric 
sideband regions (SBR), each of width 2a m and centered 
at ±5<7 m from the fitted mass, are used to represent 
the lifetime distribution of background in the SR. The 
SR and SBR are shown in Fig. 2] 

For the lifetime analysis we use the reduced proper 
time, t' = (L — 6<tl)/Pjc. The use of the reduced 
proper time ensures that only a small acceptance correc- 
tion to the lifetime distribution is needed The av- 
erage proper time resolution for this decay sample (51 fs 
for Run Period A and 43 fs for Run Period B) is small 
enough compared to the lifetime to use a binned likeli- 
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FIG. 1: The K + K ir + invariant mass distributions for the 
(a) (?!>(1020)7r + and (b) 7f*(892)°if + decay modes. The data 
is given by the points while the line gives the fit with the 
hatched region showing the fitted background level. The ver- 
tical dotted lines gives the signal and sideband regions. 



hood method. We fit the 0(1O2O)tt and K*(892)°K decay 
modes separately. For each decay mode the Run Period 
A and B data are used in a combined fit. 

For each of the run periods A and B, the t' distribu- 
tions for the decays in the SR and SBR are binned into 
two separate histograms from 0-3 ps in 60 fs bins. The 
observed number of decays in the i th t' bin is sf for the 
SR and bf for the SBR, where A refers to Run Period 
A. The t' distribution of the SBR is used as a measure 
of the lifetime distribution of background events in the 
SR. Thus the expected number of decays (rii(A)) in the 
i th t' bin of the SR for Run Period A is given by: 



rii(A) = Sa: 



and the likelihood for Run Period A is given by 



c A = U 



m(Ay 
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NA\ 



(1) 



(2) 



Sa is the total number of signal events, Ba is the total 
number of background events in the SR, and Sa + Ba — 
T,sf. The total number of events in the SBR is = 
Tiibf and cxa is the ratio of the number of events in the 
SBR to the number of background events in the SR. The 
value of a a is obtained from the fit to the invariant mass 
distribution and is very close to 1. There is a similar 
likelihood for Run Period B and the likelihood that is 
maximized is C = Ca x £-b- The fit parameters are Ba, 
Bb, and r. 

The effects of geometrical acceptance, detector, trig- 
ger and reconstruction efficiencies, and absorption are 
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FIG. 2: The f(t') correction functions for the two decay 
modes and two run periods (a) 0(lO2O)7r Run Period B; (b) 
0(lO2O)7r Run Period A; (c) K*(892)°K Run Period B; and 
(d) ~K*i892)°K Run Period A. Deviation from a flat line 
indicates the correction from a pure exponential. 



given by the fit') correction function. The f(t') func- 
tion is determined using a detailed Monte Carlo (MC) 
simulation of the experiment where the production (us- 
ing Pythia [l^]) was tuned so that the production dis- 
tributions for data and MC matched. A full coherent 
£>+ — ► K~K + tt + decay was simulated using results from 
FOCUS [3. The /(*') distributions are shown in Fig.H 

We obtain fitted lifetimes of 507.60 ± 6.46 fs and 
506.90 ±10.60 fs for the 0(lO2O)7r and K* (892)° K decay 
modes, respectively. The fit confidence levels are 2.0% 
and 0.13%, respectively. As discussed below, these low 
values are expected for the fitting technique used. The 
lifetime distribution of all decays in the SR are shown in 
Fig. |3 together with the fit and the level of background 
contained in the SR. 

Detailed studies were performed to determine the sys- 
tematic uncertainty in the lifetime measurement. 

The fit') correction reduces the fitted lifetime by 
0.80%. We studied the uncertainty in this correction. We 
verified that the MC reproduces the data Df longitudi- 
nal and transverse momenta, the multiplicity of the pro- 
duction vertex, cos 9* , and the decay length and proper 
time resolutions. A sensitive check of the acceptance and 
efficiency part of the MC correction was done using high 
statistics Kg ir + ir~ decays. Short-lived K s decays 
were reconstructed using the same analysis methods in 
the same decay region as the decays. Since the K s 
lifetime is well known, we can determine the fit') cor- 
rection in data and compare it to that obtained in our 
MC simulation. The agreement is excellent but was lim- 
ited in sensitivity by both data and MC statistics. Using 
this sensitivity as the level of the uncertainty in the fit') 
correction, we determine a systematic uncertainty due to 
this correction of ±0.69%. Possible time dependent sys- 
tematic effects were looked for by splitting the data into 
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FIG. 3: The lifetime distribution for all decays in the data 
SR (points), and the fit (histogram). The shaded distribution 
shows the lifetime distribution of the background component 
in the SR. The data and fit are shown for Run Periods A and 
B combined in the same plot for the two decay modes (a) 
0(1O2O)tt and (b) K*i892)°K. 



different time periods and comparing the fitted lifetimes. 
We also compared the separate fitted lifetimes for decays 
originating from each of the four targets. No systematic 
uncertainties were found in these comparisons. 

Our limited knowledge of the production and decay 
of the Df could contribute to a systematic uncertainty. 
This was studied using different MC simulations where 
the production parameters and the resonance substruc- 
ture of the decay were varied over reasonable ranges. 
Production systematics were also studied by splitting the 
data into different bins of total and transverse Df mo- 
menta, primary vertex multiplicity, and by comparing 
the fitted lifetimes for particles and anti-particles. We 
found no evidence of systematic effects, however, vari- 
ations of the fitted lifetimes in these studies led us to 
assign a systematic uncertainty of ±0.23% due to our 
limited knowledge of Dj production and decay. 

The systematic uncertainty due to absorption of the 
Dj and daughter particles was determined by varying 
the .D+ interaction cross-section [23| by 100% and the 
daughter particle interaction cross-sections by 50% in the 
MC. ft was also studied by comparing the lifetimes of 
decays occurring inside and outside of the target, and 
by comparing the lifetimes for decays where the Df was 
produced in the upstream half of each target with those 
produced in the downstream half of the same target. We 
determined a systematic uncertainty of ±0.38% due to 
absorption. 

In order to use the reduced proper time we must be 
able to correctly model our proper time resolution. This 
was verified by comparing the distributions for data and 
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MC and by studying splits of the data sample that can be 
sensitive to resolution effects. This included splits in t' 
resolution and reconstructed invariant mass and varying 
the t' bin width and changing the fitted range for t'. 

The lifetime distribution of the background in the SR 
is assumed to be well represented by the SBR. The un- 
certainties that arise because of this assumption were de- 
termined by a number of studies. 

The contamination from D + — * K~Tt + Tr + , and A+ — > 
K~pir + decays misidentified as K~K + n + decays were 
determined in our sample. As a fraction of the total 
D+ -> 0(1O2O)tt (L>+ -> K*(892fK) signal in the SR, 
we found the above two backgrounds contribute 0.7% 
(3.0%) and 0.2% (2.5%), respectively. The small con- 
tribution of these reflection backgrounds mean they give 
rise to insignificant uncertainties. This was verified in a 
test by using stronger particle identification and by ex- 
plicitly eliminating them by cutting out the appropriate 
mass regions. 

The background lifetime uncertainty was further inves- 
tigated by using symmetric sidebands of different widths 
(2-4cr m ), and located at different separations from the 
signal region (±4 to ±6<r m ). Variations in particle iden- 
tification and vertexing selection were used to study the 
lifetime with changes in the signal/background ratio and 
this showed no significant effect. 

Finally, an independent analysis which did not rely on 
knowledge of the background lifetime distribution was 
performed. In this analysis the data were split into fifteen 
150 fs wide reduced proper time bins from 0-2.25 ps and 
one 750 fs wide bin from 2.25-3 ps. The numbers of 
D+ -> 0(1O2O)tt+ and K*(892)°K + decays in each bin 
were determined in a mass fit and the yields for the two 
decay modes separately fitted to an exponential decay 
distribution modified by a /(£') correction function. No 
significant systematic effects were found. To account for 
the variations of the lifetime in the various background 
studies we assign a background systematic uncertainty of 
±0.16% (±1.09%) for the 0(1O2O)tt (K*(892)°K) decay 
mode. 

A mini-Monte Carlo test of the fit method was done 
to check the accuracy of the fit errors and to look for 
possible biases in the fit. We generated 1000 indepen- 
dent data samples from the best fit to the data. Each of 
these data samples were fit as real data. Distributions 
of fitted values showed a small bias in the fitted lifetime 
which is accounted for in a Fit Method systematic uncer- 
tainty (see Table HJ. This study also showed that the fit 
error underestimates the true statistical uncertainty in 
the </>(1020)tt (K*(892) K) mode by 4.8% (14.7%). The 
statistical errors on the lifetimes quoted in this paper 
have been corrected for this effect. These fit properties 
are due to neglecting statistical variations in the shape 
of the background in the SR, which is fixed in the fit 
to be equal to the measured background in the SBR as 
given by b i ' . This also leads to the small reported fit 



confidence levels. When bin-to-bin fluctuations in the 
observed b i ' are accounted for in the calculation of the 

TABLE I: Contributions to the systematic uncertainty. 



Contribution 


Systematic Uncertainty (%) 




0(1O2O)tt 


K (892)°A' 


Combined 


Acceptance 


±0.69 


±0.69 


±0.69 


Production 


±0.23 


±0.23 


±0.23 


Absorption 


±0.38 


±0.38 


±0.38 


Backgrounds 


±0.16 


±1.09 


±0.41 


Fit Method 


±0.12 


±1.07 


±0.38 


Time Scale [15J 


±0.11 


±0.11 


±0.11 


Total 


±0.85 


±1.74 


±1.00 



TABLE II: Comparison of Df lifetime measurements and 
t(D+)/t(D°) ratios. The PDG value of the D° lifetime 
is used in the ratios, except for the FOCUS ratio for which 
the FOCUS measurement is used [r^|. 



Experiment 


t(D+) fs 


r(Df)/r(D°) 


E687 [1£] 


475 ± 20 ± 7 


1.158 ±0.052 


E791 [19] 


518 ±14 ±7 


1.262 ±0.038 


CLEO II.5 [20J 


486.3 ± 15.0±g? 


1.185 ±0.039 


SELEX [21] 


472.5 ± 17.2 ±6.6 


1.152 ± 0.045 


FOCUS (this result) 


507.4 ±5.5 ±5.1 


1.239 ± 0.017 



X 2 , the fit confidence levels for 0(1O2O)tt and K*(892)°K 
become 8.8% and 16%, respectively. The goodness of fit 
is further verified by using background lifetime distribu- 
tions taken from fits to the data sideband distributions 
instead of the actual (bf' B ) entries. This fit method gives 
consistent lifetimes and fit confidence levels of 9.2% and 
18% for 0(1O2O)tt and K*(892)°K, respectively. 

We have measured the Df lifetime to be 507.4 ± 
5.5 (stat.) ± 5.1 (syst.) fs using 8961 ± 105 D+ -> 
(/)(1020)7r+ and 4680 ± 90 D+ -> K*(892)°K + decays 
from the Fermilab FOCUS photoproduction experiment. 
We also determined the ratio of our D+ lifetime to our 
D° lifetime [Hi to be 1.239 ±0.014 (stat.) ±0.009 (syst.), 
where the systematic uncertainty in the ratio is reduced 
by eliminating common systematic contributions. Ta- 
ble [n] compares our measurement with recent published 
results. This measurement represents a significant im- 
provement in precision. 
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